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Анотація: Розрахункові схеми, в яких умови закріплення можуть бути промодельовані як 
циліндричні шарніри є одні з найпоширеніших у машино та авіабудуванні. Однак точні методи 
розрахунку для таких схем при заданих умовах контактного типу відсутні. Через це створення 
алгоритму для вирішення таких задач є актуальним. Дана робота присвячена моделюванню 
нескінченного шару, що спирається на два циліндричних врізаних шарніри які паралельні 
горизонтальним поверхням шару та один одному. Шарніри представлені як циліндричні 
порожнини із заданими умовами контактного типу (нормальні переміщення та дотичні 
напруження уздовж поверхні порожнини дорівнюють нулю). На поверхнях шару напруження 
вважаються відомими. Застосовано два типу систем координат: декартова для шару та 
циліндрична - для порожнин. Граничні умови застосовуються до рівняння Ламе. За 
допомогою узагальненого методу Фур'є задача зведена до розв'язання системи нескінчених 
лінійних алгебраїчних рівнянь, яка вирішена методом редукції. Чисельне дослідження 
проведене з точністю виконання граничних умов 107 для значень напружень від 0 до 1 при 
четвертому порядку системи рівнянь. Аналіз напруженого стану проведений для опорних 
частин шару та верхньої межі шару. Виконано порівняння з результатами дослідження для 
аналогічної задачи з жорстким закріпленням опор. Відмічено, що окружні напруження в шарі 
уздовж поверхні, що контактує з циліндричним шарніром, при закріпленні контактного типу 
у порівнянні з жорстким закріпленням має протилежний знак, що може бути суттевим для 
матеріалів з різними характеристиками на розтяг-стискання. Отриманий метод може бути 
використаний для аналізу напружено-деформованого стану болтового чи заклепочного 
з'єднання в якому під час експлуатації змінився момент затягання і з'єднання з жорсткого 
перетворилось на шарнірне. 

Ключові слова: шар з циліндричними порожнинами, узагальнений метод Фур'є, умови 
контактного типу 
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1. Вступ 


Задача визначення напружено-деформованого стану є однією з найбільш важливих при 
проектуванні машин та механізмів в різних галузях авіа та машинобудування. При цьому тип 
опор у вигляді циліндричних шарнірів є найбільш розповсюдженим. 

Найбільш раціональним підходом до моделювання циліндричного шарніру є 
представлення його у вигляді порожнини або включення з заданими умовами контактного 
типу. 

Літературний пошук показує велику кількість робіт де розглядається шар з циліндричною 
порожниною або включенням які перпендикулярні до поверхонь шару [1 — 5]. Так у роботі [1] 
побудовано хвильове поле нескінченного пружного шару, ослабленого циліндричною 
порожниною. Умови закріплення — ідеальний контакт уздовж верхньої та нижньої поверхонь 
шару. Зовнішнє навантаження - розтягуюча сила прикладена до поверхні циліндричної 
порожнини в початковий момент часу. Робота [2] має аналогічну постановку задачі, з 
додатковою умовою жорсткого закріплення нижньої поверхні шару. Підхід у розв'язанні 
задачі в роботах [1, 2] також аналогічний - за допомогою інтегральних перетворень Лапласа 
та інтегрального методу Фур'є. Ці методи застосовано до граничних умов і осесиметричних 
рівнянь руху, які створюють одновимірну векторну неоднорідну крайову задачу. Недоліком 
методів, використаних в роботах [1, 2], є неможливість ефективного вирішення задач з 
декількома граничними поверхнями. Таким чином, вони можуть бути застосовані лише до 
вирішення задач дифракції хвиль. 

Розробці та застосуванню загального аналітичного рішення для функціонально 
градуйованих (ЕС) композитних ламінованих пластин із круглими вирізами за різних умов 
навантаження в площині присвячена робота [3]. Надійність 1 точність запропонованого 
аналітичного рішення встановлюється шляхом порівняння його результатів з аналітичними 
результатами, отриманими за допомогою методу комплексної змінної та методом скінчених 
елементів. 

Робота |4| присвячена знаходженню оптимального розподілу термічних напружень у 
симетричних композитних пластинах з некруглими отворами під дією рівномірного теплового 
потоку. Знайдено аналітичне рішення для ламінованих пластин з композиту з круглими 
отворами, яке потім за допомогою функції відображення розповсюджене на пластину з 
некруглими отворами. Для аналізу напружень були застосовані генетичний алгоритм (СА), 
алгоритм гравітаційного пошуку (GSA) i алгоритм Bat (BA). На жаль вищеназвані методи 
методи також є наближеними що не гарантує точність кінцевого результату. 

В роботі [5] за допомогою інтегральних перетворень Вебера-Орра розглянуто кручення 
пружного півпростору з вертикальною циліндричною порожниною коаксіальним штампом, 
що обертається під дією крутного моменту навколо осі. В роботі запропоновано два нові 
методи зведення задачі до регулярних інтегральних рівнянь другого роду, що дозволило 
отримати майже точний розв'язок задачі для всіх значень параметрів. 

Загальним недоліком усіх вищезгаданих робот [1 — 5] є Te, що вони розроблені для 
вирішення задач з поперечним розташуванням порожнини або включення, і без дуже суттєвого 
доопрацювання їх не можна використати для розв'язання задачі для шару з поздовжніми 
циліндричними порожнинами. 

Інший підхід до розв'язання задач для шару, що має повздовжнє контактне спирання це 
застосування чисельних методів та методів будівельної механіки. Аналіз існуючих робіт 
показав, що в більшості робот використовувався метод скінчених елементів [6, 7). Прикладом 
таких робіт є дослідження [8] де за допомогою методу скінчених елементів проведено аналіз 
напружено-деформованого стану півпростору, підкріпленого плитою. Значна відмінність в 
постановці задачі даного дослідження є вертикальне розташування армованої оболонки 
циліндричної порожнини. Недоліком даного класу методів є те, що при використанні методів 
будівельної механіки вони допускають суттєве змінення моделі, метод скінченних елементів є 
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наближеним, не враховує нескінчені границі тіла та суттєво залежить від точності 1 розмірів 
побудови сітки. Тому дані методи не гарантують точність та адекватність результатів [9]. 

Застосування методів, що використовуються для розрахунку композитних матеріалів [10 
- 13], дозволяє врахувати нелінійність моделі. Циліндричні шарнірні опори при цьому 
можливо моделювати у вигляді абсолютно жорстких елементів. Дослідження [10] присвячено 
визначенню напружено-деформованого стану під дією поздовжньої динамічної сили, 
прикладеної до одного з двох з'єднаних внахльост стрижнів різної довжини. Для моделювання 
з'єднання використовувалась модель клейового з'єднання Голанда-Рейснера. Напружено 
деформований стан при поперечному ударі розглянуто в роботі (11). Предметом дослідження 
при цьому виступали багатошарові конструкції до яких прикладено динамічне навантаження. 
Рішення було отримано за допомогою теорії двовимірної дискретної структури, при цьому 
вектор переміщень кожного шару був представлений у вигляді степеневого ряду. Крім 
теоретичних проводилися експериментальні дослідження, які підтвердили коректність 
результатів. Також дослідженню ударного навантаження присвячена робота [12], в якій 
створено метод оцінки міцності багатошарового склопакету при зіткненні з птахом. Крім того, 
в роботі оцінювалася величина надлишкового тиску під час ударного навантаження на 
багатошаровий склопакет. В роботі [13] проведено дослідження термонапруженого стану 
багатошарового склопакету. При цьому термічне навантаження промодельоване у вигляді 
міжшарових плівкових джерел тепла, а склопакет розглянуто як незамкнута циліндрична 
багатошарова оболонка постійної товщини. 

Суттєвим недоліком даного класу методів є неможливість врахування циліндричних опор 
(нульових переміщень в зонах контакту з циліндричними шарнірами). 

Застосування аналтико-числового узагальненого методу Фур'є [14] є найбільш 
ефективним для врахування поздовжніх неоднорідностей в шарі. Цей метод дозволяє 
отримувати розв'язок рівняння Ламе для системи тіл, в якій використовується декілька 
локальних систем координат різних типів одночасно. Використовуючи даний метод отримано 
напружено - деформований стан в пружному циліндрі з циліндричними порожнинами | 15, 16] 
або циліндричними включеннями [17], а також для півпростору зі сфероїдальною порожниною 
[18]. При цьому використовувались локальні циліндричні системи координат та формули 
переходу базисних рішень між ними. Для шару зі сфероїдальною порожниною 
використовувалась декартова система координат та сфероїдальна. 

Формули переходу між декартовою та циліндричною системою координат застосовані 
для півпростору з циліндричною порожниною в роботі [19], для шару з однією циліндричною 
порожниною в роботах |20, 21], для шару з суцільним циліндричним включенням в роботі [22], 
а також для шару з циліндричною товстостінною трубою в роботі [23]. 

Загальним недоліком робіт [19-23] є те, що в них не застосовуються формули переходу 
між локальними циліндричними системами координат. Даний недолік не дозволяє 
розв'язувати задачі з декількома циліндричними неоднорідностями. 

В подальшому розвитку робіт [19—23] була розв'язана задача для шару з двома суцільними 
циліндричними включеннями [24], але в даному випадку включення не можуть розглядатися 
як опори. В роботі |25| розв'язана задача визначення напружено деформованого стану в шарі 
з двома циліндричними врізаними та жорстко поєднаними з шаром опорами. Але в цій роботі 
на поверхнях порожнин застосовується граничні умови у вигляді переміщень, що не дозволяє 
використовувати їх для вирішення задач з умовами контактного типу. 

В той же час при створенні розрахункових схем машинобудівних та авіаційних 
конструкцій часто застосовується схема, де шар спирається на два та більше циліндричних 
шарніри. Розв'язок задачі знаходження напружень за допомогою аналітико-чисельних 
високоточних методів (наприклад узагальненого методу Фур'є), як показана вище, не 
розглядався раніше. 

Таким чином, задача розв'язку задач контактного типу для шару, що спирається на два 
циліндричних шарніри за допомогою узагальненого методу Фур'є є актуальною. 
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2. Постановка задачі 


Пружний однорідний шар розташований на двох циліндричних врізаних в нього шарнірах. 
Ці шарніри було представлено у розрахунковій схемі як порожнини із заданими на них 
умовами контактного типу (рис.1). Радіуси порожнин позначені Rp, де р - номер порожнини. 


o z) 


Рис. 1. Шар з двома врізаними циліндричними опорами. 


Локальна система координат для шара - декартова (х, у, <), для порожнин - циліндрична 


(Pp, Фр» 2). Відстань до верхньої та нижньої меж шару відповідно y=htay=—h. 
При вирішенні задачі знайдено розв'язок рівняння Ламе: 
Au + (1 - 20)71Vdivu = 0. 
На верхній Ta нижній межах шару задано напруження, відповідно F U(x, 2)у=- = 
20 > 20 
F, (х, 2), FU(x, 2) „= = F; Œz), 
де И — переміщення в шарі; 
> © > oR, ðR 1f> > Р 
ЕЦП = 2:6: [i аїр + on U+ 5 (п x rotU)| - оператор напруження; 


= h) > h) > h) > 
F? (x, z) = 1098, + os Ne; + тр, (1) 
2 h) > h) > h) -> 
F?(x,z) = a, + al г, + та, (2) 


- відомі функції. 
На поверхнях порожнин р = 1, p=2 задані нормальні переміщення та дотичні 
напруження 


А = U (pa). 
Tpplop=Rp — aO) (3) 


Toz|pp=Rp 7 г” (Ф, 2) 


Всі задані функції прийнято вважати швидко спадаючими від початку координат по осі Z 
и х для шару та осі Z для порожнин. 
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3. Метод розв'язання 


Базисні розв'язки рівняння Ламе для декартових 1 циліндричних систем координат було 
взято у вигляді описаному в роботі | 141: 


1. (GY, Z; A, u) = МФ eiGztuxtry, 
Втр, 0,2; A) = NP In (Apeiaz+me), 
Зкт р, $, 220 = NP (відп dy" Kn (А) ‚ант k = 1,2,3; а) 
МӨ =17; МӨ = DEP +17(у ); М9 = о); 


кезе Диане о PE 
МР = srot( és” ); у= 4/42 +u, – 0 < 24, < ә, 


де Їн (2), Km(x) - модифіковані функції Бесселя; Rkm» 5 т — Внутрішні Й зовнішні 
розв'язки рівняння Ламе для циліндру відповідно; С ), Шо в розв'язки рівняння Ламе для 


шару; с - коефіцієнт Пуассона. 
Розв'язання задачі представлено у вигляді запропонованому в роботі [25]: 


г-у У | 2 Bem A) кт (Рр» Фр» ZA) dA + 


p=1 k= 
+ Diaa Го Го (HA Обь ув Aw) + А, р) HE оу, 25 A, u) ) анал, (5) 


де т (рь, Фр» 27), Ох, у, Zz; À, и) і (x, у, Z; À, и) - базисні розв'язки, задані формулами 


(4), а невідомі функції Н, (А, и), Ну (А, и), Be т (2) знайдені з граничних умов (1) - (3). 

Перехід між базисними розв'язками в різних системах координат було здійснено за 
допомогою формул отриманих в роботі [14]: 

- ДЛЯ во від базисних розв'язків Sin циліндричної системи координат до 


розв'язків шару 10 2 (при y>0) 1 йо (при у<0) 


-і)т n р а 
Зкт( р Pp» z; 1) = ~ > Г OT . егінхріуУр . 10) = k= 
Sam(Pp бша a) 7-72 > | we |1 Ет р — ae 2° Yp F дай) + 


ет |ИХрЯУУр du 


£4u(1 — ous?) 2—6, (6) 
деу = 4/2 + и?, ох (А, и) = E, m = 0, +1, +2,...; 


1, 3; 


- для переходу від базисних розв'язків 0001 С ) шару до розв'язків Rew циліндричної 


системи координат 


o0 


ui (x,y,z) = е ИХрРЯУУр З > (i . w)” Rem, (k = 1, 3); 


тезою 
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a (x, y, Z) — ер-УУр . 2 la . wz)” . 472 ((т . и + Yp О 22) . Rim + 
те-ою 
+y Rom + 4u(1—0)Ram)] , a 
де 


> 


= Drm( Pp, 4) - є (Фотій); By (ру) = ëp (Ар) + i- р). (2% + ê); 


з А Я МР 4(0 - 1) 
bon(p, 1) = ер і (40 Е: 3) I n(Ap) + АРІ п (Яр) + еф! 7 хі np ) е Ар г Ор зб 
+é,iApI'n (Ap) ; 

ваз (р, А) = — lé, Ор) + by i 1 (ар); 


en: ёф, ё, — орти в циліндричній системі координат; 
- для переходу від базисних розв'язків циліндру з номером р до розв'язків циліндру з 
номером 4 


Skm(Pp Pp Zi A) = Ух. „Бан (Pa) Ое, 2,218 


brs. (ра) = (-1)"Kin- aAa): ної б, „(ра,4); 
Бит (ра) SL) Ke (Ара). E eee: баа (во) (8) 


В (pq) = = (– "(к т- (ба) 5, nlPq A) =, 
| [Кина (бра) + Kin—n—1 A€pa) Ви (р А) ve NM Hipi 


де дра - КУТ між віссю Xp та відрізком вр, Ки (0) = (sign(x))” «Ку, (\xI). 


Для врахування граничних умов на верхній Ta нижній межах шару застосовано оператор 
напруження для правої частини функції (5). Отриманий вектор дорівнює при у = h заданому 
FP (х, 2), а при у =—h заданому Б? (х,2), які попередньо представлені через подвійний інтеграл 
Фур'є. 

Для базисних розв'язків Sk m (Pp Pp» 2; 1) застосовано формули переходу (6), представши 
їх у декартовій системі координат через базисні розв'язки WO? (x, у,2;Л,и) при y=h і 
(+ т А РУБ ад 
10 (x, у, Z; À, и) при y =—h. Після врахування граничних умов на верхній і нижній поверхнях 
шару отримано шість рівнянь (по одному на кожну проекцію) з 12-ма невідомими НІ (4, и), 

2 
А. (2, и), Bem (А). 

Із цієї системи рівнянь знайдено Н, (A, и) i Ну (A, и) через Be mA): 

Для врахування граничних умов контактного типу на циліндрах, праву частину (5), за 
допомогою формул переходу (6) та (7), було представлено в локальній циліндричній системі 
координат кожної порожнини через базисні розв'язки Кі т» Sk m- Після цього до цих базисних 
розв'язків, при = 2 та К = 3, було застосовано оператор напруження. Отриманий вектор, при 
Pp=Rp, прирівняно заданим (3), представленим через ряд та інтеграл Фур'є. Таким чином, 
отримано додатково ще шість рівнянь які задовольняють граничним умовам на поверхні 
циліндричних шарнірів. 

У результаті для кожного циліндра отримано по три нескінчених системи лінійних 


інтегро-алгебраїчних рівняння відносно BY т (A), які містять Нк(Я, и) i Нк (А, и). Після 
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виключення з цих рівнянь знайдених раніше Ну (Я, и) і Ну (Л, и) через BY (A) та звільнення їх 
від рядів по т i інтегралів по A, отримано 6 нескінчених лінійних алгебраїчних рівнянь другого 


роду для визначення невідомих go mA). 


Знайдені невідомі В, т (2) підставлено у вираз для Нк (Я, и) i A kà и). 

Таким чином під час розв'язку задачі було знайдено всі невідомі виразу (5). 

До отриманих нескінченних систем рівнянь застосований метод редукції. Точність 
виконання граничних умов при чисельному дослідженні показали високу збіжність 
розв'язання цієї системи лінійних алгебраїчних рівнянь. 


4. Чисельні дослідження напруженого стану 


У пружному ізотропному шарі (рис.1) розташовано дві циліндричні порожнини. 
Коефіцієнт Пуассона шару (ABS пластик) с = 0,38, модуль пружності E= 1700 Н/мм2. 
Геометричні параметри моделі: Ro= К; = 5 мм, h = 12 мм, h= 12 мм, 01. = 0, а1з= п. Відстань 
до порожнин /.2 = L13 = 30 мм. 

На верхній межі шару задано нормальні напруження у вигляді одиничної хвилі а (x,z) = 
—108 - (22 + 102)-? - (x? + 102)-? і нульові дотичні напруження б ) = = = 0, на нижній 
межі шару задано напруження, що дорівнюють нулю: с х, 2) = т х, 2) = т) (х,2) = 0. 

Нескінчена система була зрізана по параметру m=4 (кількість членів ряду Фур'є і порядок 
системи рівнянь). 

Обчислення інтегралів виконано квадратурними формулами Філона. Точність виконання 
граничних умов при зазначених т і заданих геометричних параметрах не менше ніж 10“ при 
значеннях від 0 до 1. Це відповідає роботі [20], в якій проведений грунтовний аналіз збіжності 
ров'язків відносно різних значень M для шару з порожнинами. 

Отримані результати наведені на рисунках 2- 5. Для порівняння та аналізу результатів на 
рисунки також додано графіки які було надані в роботі [25] для аналогічної задачі з жорстким 
закріпленням опор. Ці лінії позначені цифрою 1. Усі геометричні розміри, параметри матеріалу 
та завнішні навантаження, для коректності порівнянь, взяті однаковими. 

На рис. 2 наведено графік розподілення нормальних напружень уздовж верхньої та 
нижньої меж шару при 270, де оу - задані напруження, а Oy напруження, що виникають під 
дією зовнішнього навантаження. 


ЕВА 
ПО «ОП Е 


-40 -30 -20 -10 0 10 20 30 40 


Напруження, МПа 


Відстань, мм 


Рис. 2. Напруження на верхній і нижній поверхнях шару: | - задані oy; 2 - Oy на 
верхній межі при жорсткому закріпленні опор; 3 - о, на верхній межі при закріпленні 
контактного типу; 4 – Oy на нижній межі при жорсткому закріпленні опор; 5 - 0, на нижній 
межі при закріпленні контактного типу. 
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Як видно з графіків: зміна умов закріплення не впливає на характер розподілення напружень. 
Величини максимальних значень на верхній поверхні шару (лінії 2 та 3 на рис. 2) відрізняються 
приблизно на 13%. 

На рис. 3 зображено напружень Op, Og та с, на поверхні опори, розташованої праворуч 
(р = 1) при 2-0. 

Як видно з графіків, характер розподілення напружень о, та с, для жорсткого закріплення 
та закріплення контактного типу також аналогічний. В той же час значання максимальних 
величин напружень для шару, який розташований на опорах з закріпленням контактного типу 
більші більш ніж вдвічі. 
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хо 
o X 


> 
"> 


0 л/4  т/2 31/4 т 51/4 61/4 ТА 27 
Кут повороту, радіан 
умови контактного типу 
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Nn 
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о 
D 
N 


— — жорстке з'єднання умови контактного типу 


в 
а – напруження Op; 6 - напруженняоф; в - напруженняо, 
Рис. 3. Напруження на поверхні порожнини р = 1. 
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В той же час напруження Og при закріпленні контактного типу у порівнянні з жорстким 
закріпленням має протилежний знак (дивись рис.36), що може бути суттевим для матеріалів з 
різними характеристиками на розтяг-стискання, так як що дана модель може бути використана 
як модель болтового чи заклепочного з'єднання в якому під час експлуатації змінився момент 
затягання 1 з'єднання з жорсткого перетворилось на шарнірне. 


5. Висновки 


В роботі розв'язана задача теорії пружності для шару розташованого на двох кругових 
циліндричних врізаних в нього опорах. Опори представлені у вигляді циліндричних порожнин 
з умовами контактного типу на їх поверхнях. Задачу розв'язано за допомогою аналітико- 
числового узагальненого методу Фур'є. При цьому розв'язок виконаний за допомогою 
зведення задачі до нескінченної системи лінійних алгебраїчних рівнянь, яка дозволяє 
застосувати до неї методу редукції. Чисельні дослідження показали високу точність виконання 
граничних умов, що слугує додатковим підтвердженням точності розрахуннків. 

Представлений порівняльний аналіз 3 результатами дослідження в якому 
використовувались інші умови закріпленння показує, що зміна умов з жорсткого закріплення 
на закріплення контактного типу значним чином впливає на величину 1 знак напруженого 
стану тіла. 

В роботі представлено метод, що дозволяє отримати напружено-деформований стан для 
нескінченного шару, що спирається на два циліндричних врізаних шарніри які паралельні 
горизонтальним поверхням шару та один одному. В подальших дослідженнях планується 
модифікація даного методу для збільшення кількості опор та врахування можливих циліндричні 
неоднорідностей між опорами. 
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Analysis of the stress state of a layer supported by two cylindrical swivel joints 
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Abstract: Free body diagrams in which connections can be modeled as cylindrical swivel joints 
are among the most common in machine and aircraft engineering. However, there are no exact 
calculation methods for such schemes under given contact type conditions. Because of this, the 
creation of an algorithm for solving such problems is actual goal. The current investigation is devoted 
to the modeling of an infinite layer supported by two cylindrical swivel joints with parallel horizontal 
layer surfaces and each other. The swivel joints are represented as cylindrical cavities with specified 
contact type conditions (normal displacements and tangential stresses along the surface of the cavity 
are equal to zero). The stresses on the surfaces of the layer are assumed to be known. Two types of 
coordinate systems are used: Cartesian for the layer and cylindrical for cavities. Boundary conditions 
are applied to the Lamé equation. Using the generalized Fourier method, the problem is reduced to 
solving a system of infinite linear algebraic equations, which is solved by the reduction method. The 
numerical study was carried out with the accuracy of fulfilling the boundary conditions 10-3 for stress 
values from 0 to 1 at the fourth order of the system of equations. The analysis of the stress state was 
carried out for the supporting parts of the layer and the upper boundary of the layer. A comparison 
was made with the results of the study for a similar task with fixed supports. It was noted that the 
circumferential stresses in the layer along the surface in contact with the swivel joint, when the contact 
type is fastened, compared to the rigid fastening, has the opposite sign, which can be essential for 
materials with different tensile-compressive characteristics. The current method can be used to 
analyze the stress-strain state of a bolted or riveted connection in which the tightening moment 
changed during operation and the connection turned from rigid to hinged. 

Keywords: layer with cylindrical cavities, generalized Fourier method, contact-type conditions. 


